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Abstract. Astrophysical space missions deliver invaluable information about our universe, stellar dynamics of our galaxy, 
and motion of celestial bodies in the solar system. Astrometric space missions SIM and Gaia will determine distances to stars 
and cosmological objects as well as their physical characteristics and positions on the celestial sphere with microarcsecond 
precision. These and other space missions dedicated to exploration of the solar system are invaluable for experimental 
testing of general relativity. Permanently growing accuracy of space and ground-based astronomical observations require 
corresponding development of relativistic theory of reference frames, motion of celestial bodies, and propagation of light/radio 
signals from a source of light/radio to observer. Such theory must be based on Einstein's general relativity and account for 
various relativistic effects both in the solar system and outside of its boundary. We describe a hierarchy of the relativistic frames 
adopted by the International Astronomical Union in 2000, and outline directions for its theoretical and practical extentions 
by matching the IAU 2000 reference frames in the solar system to the cosmological Friedman-Robertson- Walker reference 
frame and to the frames used in the parametrized post-Newtonian formalism. 
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1. INTRODUCTION 

The role of inertial reference frames in astronomy has been recognized early by both theorists and observers. As- 
trometric and navigation data rely on observations done in a local inertial frame and transformed to a global inertial 
frame. The choice of the local reference frame is usually driven by instrumental considerations while construction of 
the global frame is based upon deeper theoretical grounds. This paper deals with the latter subject. 

The original approaches to construct reference frames in astronomy were completely based on the concepts of 
the Newtonian gravity and Euclidean absolute space and time [1-3]. Modern astrometry and navigation, however, 
are operating at the angular resolution already exceeding 1 milliarcsecond [4], and in some cases the resolution 
has achieved the level of 10 microarcseconds (/ias) [5]. At this level the primary gravitational theory must be the 
Einstein's general relativity (GR) with a corresponding replacement of the Euclidean space and time by the four- 
dimensional Riemannian space-time manifold endowed with the null-cone structure. In other words, the theoretical 
basis of modern astrometry and navigation must be entirely relativistic [1-3, 6, 7]. Recognition of this fact is rapidly 
spreading in the astrometric and navigation community especially after the successful completion of the HIPPARCOS 
mission [8], materialization of the International Coordinate Reference Frame in the sky based on radio quasars [9, 10], 
development of technologically new ideas and projects both in space (GP-B [11], FAME [12], DIVA [13], SIM [14], 
Gaia [15], VSOP [16]) and in ground-based (IERS [17], VLBA [18], SKA [19]) astrometry, the discussion [20] and 
the adoption [21] by the XXIV th General Assembly (GA) of the International Astronomical Union (IAU) of a new 
post-Newtonian framework for the reference frames in the Solar system. New motivations for extention of the IAU 
framework appeared after original discovery in 1998 [22] and recent, more robust experimental confirmation [23] of 
the, so-called, Pioneer anomaly observed in the hyperbolic orbital motion of spacecrafts (especially Pioneer 10 and 11) 
in deep space of the solar system. This anomaly may reflect a hidden influence of the gravitational field of our Milky 
Way galaxy and/or cosmological Hubble expansion on the local dynamics inside the solar system [24]. Furthermore, a 
microarcsecond astrometry can deliver a unique opportunity to measure the Galaxy and cosmological parameters (and 
perhaps gravitational waves from the early universe [25, 26]) from angular measurments of the position and proper 
motion of quasars in the sky. Thus, the global reference frame used for interpretation of astrometric and/or navigation 
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data should properly incorporate the galactic [27] and cosmological effects [28-30]. 

In what follows we shall describe in more detail the initial theoretical motivations that stimulated the development 
of modern relativistic framework for the reference frames in the solar system and new ideas driving the development of 
the theory beyond 2000. We also provide the reader with a brief history of the development of the relativistic reference 
frames (section 2). Main notations are explained in section 3. Then, in section 4, we give a description of the IAU 1991 
reference frame framework and compare it in section 5 with the present-day formulation adopted by the XXIVth GA of 
the IAU 2000. Matching the IAU 2000 framework with the parametrized post-Newtonian formalsim being extensively 
used in experimental gravity research, is given in section 6. Finally, we match the reference frame of the expanding 
cosmological model to the IAU 2000 resoltuions in section 7. 

Transition from the Newtonian concepts towards more profound theoretical relativistic approach for the construction 
reference frames in modern astrometry and celestial mechanics can be traced back to the period 1975-1992 and is 
associated with discovery of a binary pulsar PSR 1913+16 by Hulse and Taylor [31] and the breakthrough in the 
solution of two major theoretical problems: (1) development of a self-consistent framework for derivation of the post- 
Newtonian equations of translational and rotational motion of self-gravitating extended bodies in the solar system [32- 
35], and (2) development of the higher-order post-Newtonian celestial mechanics of binary pulsars [36-38]. The main 
concern of theorists working on relativistic celestial mechanics of the solar system bodies was related to the problem of 
unambiguous interpretation of astronomical measurements and separation of small coordinate (gauge) perturbations 
from the real physical relativistic effects. Later on, the question regarding the best choice of the gauge conditions 
imposed on the metric tensor components arose, for a number of gauges were used in calculations by various groups, 
so that there were some arguments about their advantages and disadvantages as well as about the uniqueness of the 
results [39] . But, probably the most serious was the problem of construction of the local geocentric (and planetocentric) 
reference frame in the post-Newtonian approximation(s). Apparently there existed principal difficulties in solving 
these problems. This was because the simplest post-Newtonian approximation of the Earth's gravitational field by 
the Schwarzschild spherically-symmetric solution could not be considered as accurate enough due to the noticeable 
perturbations to the metric tensor from the Earth's rotation and oblateness as well as from the tidal forces produced 
by the Moon, Sun, and other planets. It was also recognized [33] that the well-known procedure of construction of the 
Fermi normal coordinates in the vicinity of a geodesic world-line of a test particle [40, 41] implemented by Fukushima 
[42], and Ashby and Bertotti [32] for construction of the geocentric coordinate system can not be applied due to: (1) 
the ambiguity in choosing the background space-time manifold, and (2) deviation of the Earth's center-of-mass world- 
line from the geodesic motion because of coupling Earth's internal multipole moments with tidal gravitational field of 
external bodies. 

Lack of a firm theoretical approach to the problem of the construction of a local inertial frame in the neighborhood of 
a massive, extended celestial body led theorists (see, e.g., [43^-8] and references therein) to use a simple (Euclidean- 
like) spatial translation of the origin of the barycentric reference frame to the Earth's geocenter in order to make 
a geocentric reference frame at the post-Newtonain approximation. However, though such a procedure is allowed 
in general relativity because of the coordinate freedom, the Euclidean translation to the geocenter produces large 
relativistic effects having pure coordinate origin that are unobservable [49]. In addition, in such frame the geometric 
shape of the Earth moving along the elliptic orbit undergoes the Lorentz and gravitational contractions which must 
be compensated by spurious internal stresses in Earth's interior in order to prevent the appearance of unphysical 
deformations of its crust. Similar problems were also met in developing general relativistic celestial mechanics of 
binary pulsars [38, 50]. 

Recent technological developments in drag-free satellites, clocks, lasers, optical and radio interferometry [51, 52] 
make it necessary to elaborate the Will-Nordtvedt parametrized post-Newtonian formalism [43] in order to incorporate 
relativistic effects in local frames constructed around self-gravitating and extended bodies [53-55]. The domain of 
applicability of the IAU relativistic theory of reference frames [21] should be also extended outside the boundaries of 
the solar system [28]. 

New generation of astrometric satellites which include SIM [14], and Gaia [15](recently adopted as a cornerstone 
mission of ES A) requires an absolutely new approach for unambiguous interpretation of astrometric data obtained from 
the on-board optical interferometers. SIM and Gaia complement one another - SIM excels at studying bright and faint 
sources to high accuracy, while Gaia will obtain astrometric measurements of a large number of stars. SIM represents 
an entirely new measurement capability in astronomy and a revolution in astrometry that will exceed all previous 
measurements by a factor of 100-1,000. For faint stars, SIM's capability will exceed that of upcoming missions such 
as Gaia by a factor of 25-75 or more, and will aproach to 1 /las. At this level of accuracy the problem of propagation of 
light rays must be solved with taking into account relativistic effects generated by non-stationarity of the gravitational 
field of the solar and other astronomical systems emitting gravitational waves [28]. That will allow us to integrate the 



equations of light propagation in gravitational field of light-ray-deflecting bodies without artificial assumptions about 
position and velocity of the bodies from which it deflects light. Astrometric resolution in 1 fi as changes the treatment 
of the parallax, aberration, and proper motion of observed celestial objects which becomes much more involved and 
requires better theoretical definition of reference frames on curved space-time manifold [58-60]. This problem is 
highly releveant to unambigious physical interpretation of the gravitomagnetic precession of GP-B gyroscopes [11] 
which is measured relative to a binary radio star IM Pegasi [56] that is pretty close to the solar system (has parallax) 
and moves rather fast (has large proper motion). 

The other motivation for improvement of the theory of reference frames within framework of general relativity 
is associated with continuously growing precision of the Global Positioning Systsem (GPS) that is a satellite-based 
navigation system made up of a network of 24 satellites placed into orbit by the U.S. Department of Defense and 
orbiting the earth about 12,000 miles above its surface. Ephemeris data tells the GPS receiver where each GPS satellite 
should be at any time throughout the day. Each satellite transmits ephemeris data showing the orbital information for 
that satellite and for every other satellite in the system. The GPS uses accurate, stable atomic clocks in satellites and 
on the ground to provide world-wide position and time determination. The best available current precision of the 
system in range measurements is about 1 cm, and it may be improved by a factor of 10 with making use of radio links 
of multiple frequencies. This precision suggests that the GPS data processing must include a number of relativistic 
effects [61] and definition of geodetic parameters must rely upon general theory of relativity [62]. 

Additional motivation for improving relativistic theory of reference frames can be relegated to the problem of 
calculation and detection of the incoming gravitational-wave signals of the space interferometric gravitational wave 
detectors like LISA [63] which would consist of three satellites flying apart in space separated by distances of order 
1 .5 x 10 6 km in an equilateral triangle. The center of LISA'S triangle will follow Earth's orbit around the Sun, trailing 
20 degrees behind. The spacecraft will carry instruments to track each other and measure passing gravitational waves. 
Each of these instruments is made up of two optical assemblies, which contain the main optics, lasers, and a free- 
falling gravitational reference sensor. LISA will detect gravitational wave sources from all directions in the sky. These 
sources will include all the thousands of compact binary systems containing neutron stars, black holes, or white dwarfs 
in our own Galaxy, plus merging supermassive black holes in distant galaxies. However, the detection and proper 
interpretation of the gravitational-wave signals can be done only under the condition that all coordinate-dependent (or 
gauge-dependent) effects are completely understood and subtracted from the signal. This is especially important for 
observation of gravitational waves from very distant sources located at cosmological distances because the Hubble 
expansion of our universe will affect the propagation of the waves. One should also not forget that LISA would fly 
in the near-zone of the Sun which is as a source of low-frequency gravitational waves produced by oscillations of its 
interior (known as g-modes [64]). Motion of the LISA satellites and propagation of light rays along the baseline of 
LISA in the field of such oscillating g-modes must be carefully studied in order to provide a correct interpretation of 
observations as well. 

Many other relativistic experiments for testing general relativity and alternative theories of gravity also demand 
advanced theory of relativistic reference frames (see, for example, [55, 65-67] and references therein). Two space 
laser interferometric missions - ASTROD [68] and LATOR [69] - are dedicated to measure general relativity in 
higher-order, non-linear regime. 

2. BRIEF HISTORICAL REMARKS 

Historically it was de Sitter [70, 71] who worked out a relativistic approach to build reference frames in general 
relativity. He succeeded in the derivation of relativistic equations of motion of the solar system bodies and discovered 
the main post-Newtonian effects including the gravitomagnetic perihelion precession of planetary orbits due to the 
angular momentum of the Sun (spin-orbital momentum coupling). He also discovered the geodetic (de Sitter-Fokker) 
precession of the lunar orbit due to the orbital-orbital momentum coupling that has been verified with the precision of 
< 1% [72-74]. These effects are essential in the present-day definition of the dynamic and/or kinematic rotation of a 
reference frame [21, 35]. 

Later on, Lense and Thirring [75] gave a more general treatment of the gravitomagnetic dragging of a satellite 
orbiting around a massive rotating body but they hardly believed that the effect can be measured in practice. It was 
Ginzburg [76] who first realized that the Lense-Thirring precession can indeed be measured using the earth-orbiting 
artificial satellites and paid particular attention to the problem of separation of coordinate and physical effects. But only 
recently the experimental verification of the Lense-Thirring effect has come about with a geodetic satellite LAGEOS 
[65, 66, 77]. In 1959, Pugh [78] proposed testing general relativity by observing the precession of a gyroscope 



with respect to a distant star in a drag-free satellite. Schiff [79] made calculations and came to the conclusion that 
a gyroscope in polar orbit at 400 miles should turn with the Earth through an angle amounting after one year to 
42 milliarc-seconds due to the Lense-Thirring precession. An interesting generalization of Schiff 's calulations was 
provided by Teyssandier [80] who worked out a realistic geophysical model that enables us to evaluate the contribution 
of the earth's nonsphericity to the Lense-Thirring precession of a gyroscope. The gravitomagnetic frame-drag effect 
is being currently tested by Everitt's GP-P team [11]. Terrestrial-based measurment of the gravitomagnetic precession 
in a laboratory at the south pole station was proposed by Braginskii et al. [81] and this experiment may be feasible 
with the current-day technologies. Frame dragging effect caused by translational motion of a massive body has been 
recently measured in VLBI experiment [5] that observed a minute relativistic correction to the static-field Shapiro time 
delay, caused by a moving gravitational field of Jupiter [82, 83]. 

Ginzburg's paper [76] and the success of the soviet space program motivated Brumberg to develop a post-Newtonian 
Hill-Brown theory of motion of the Moon [84, 85] and major planets [86] in the solar barycentric coordinate system 
(see also [45, 87]). Baierlein [44] extended Brumberg's approach accounting for the eccentricity of Earth's orbit. A 
different approach to the problem of motion of the Moon and the construction of reference frames in general relativity 
was suggested by Mashhoon and Theiss in a series of papers (for a review see [88, 89] and references therein). Instead 
of making use of the standard post-Newtonian approximations (PNA) they developed a post-Schwarzschild treatment 
of gravitomagnetic effects in the three-body problem. This approach treats gravitational field of the Earth in terms 
of the exact Schwarzschild solution on which small tidal perturbations from external bodies (Sun, Moon, etc.) are 
superimposed. It allowed them to discover that the validity of PNA is restricted in time due to the appearance of 
time-divergent integrals, and the geodetic and gravitomagnetic precessions are parts of the more general phenomena 
involving a long-term relativistic nutation known as the Mashhoon-Theiss anomaly [90]. The same authors also 
introduced a precise definition of the local geocentric frame with the Earth considered as a massive monopole particle 
moving on a background space-time manifold [89, 91], thus, extending the concept of the Fermi normal coordinates 
usually introduced in the neighborhood of a massless observer [92]. 

A more general approach to the problem of construction of reference frames in general relativity was initiated in 
our papers [33] (see also [93]), where the matched asymptotic expansion technique and a multipolar decomposition of 
gravitational fields of bodies were employed. This approach has been further refined and used for the development of 
the, so-called, Brumberg-Kopeikin (BK) formalism for building relativistic astronomical frames in the solar system (or 
in any other N-body system). This formalism relies upon Einstein's equations, which are solved for the construction of 
local geocentric and global barycentric reference systems whose properties are determined by the specific form of the 
metric tensor. The metric-tensor matching technique is used for setting up relativistic transformations between these 
two systems [34, 35, 87] which generalize the Lorentz transformation of special relativity to the case of presence of 
gravitational field. 

Damour, Soffel and Xu (DSX) extended the Brumberg-Kopeikin theory by protracting the Newtonian definitions of 
the multipole moments of gravitating bodies used by Brumberg and Kopeikin to the post-Newtonian approximation 
of general relativity [94-96]. Combined BK-DSX formalism made impact on the formulation of resolutions of the 
General Assembly of the IAU 1991 on reference frames [20, 97]. Further lively discussion on the reference frames 
among members of a joint committee on relativity [98] of the Bureau International des Poids et Mesures (BIPM) 
and of the International Astronomical Union made the complete BK-DSX theory accepted by the XXIV-th General 
Assembly of the IAU 2000 [21] as a basic framework for setting up relativistic time scales and astronomical reference 
frames in the solar system. 

Further development of the reference frame theory goes to the extention of the BK-DSX theory to alternative theories 
of gravity [54, 55] and matching this theory with the cosmological solution which accounts for the Hubble expansion 
[29, 99, 100]. The BK-DSX theory has been supplemented with the precise mathematical solution of the problem of 
propagation of test particles and light rays in multipolar gravitational fields of moving bodies [57-60, 101-108] that 
allows us to apply that theory for the practical aims of astrometry and navigation in the solar system. 

Recognizing these goals and taking into account that relativitistic questions played an important role in the work of 
several IAU Commissions and IAU Working Groups for an extended period of time, the XXVI-th General Assembly of 
the IAU, held in Prague in Aufust 2006, has decided to establish a new IAU Commission 52 "Relativity in Fundamental 
Astronomy" (RIFA). The general scientific goals of the new commission are 

• clarify geometric and dynamical concepts of fundamental astronomy within relativistic framework; 

• provide adequate mathematical and physical formulations to be used in Fundamental Astronomy; 

• deepen understandings of the above results among astronomers and students in astronomy; 



• promote research needed to acomplish these tasks. 
This paper pays specific attention to all these issues. 



3. NOTATIONS 

Latin indices takes values 1,2,3, and the Greek ones run from to 3. Repeated indices indicate the Einstein summation 
rule. The Kroneker symbol (unit matrix) is denoted 5y = diag(l, 1,1) and the fully anti-symmetric symbol of Levi- 
Civita Eijk is defined in such a way that £123 = 1. The Minkowski metric is r] a p = diag(— 1, 1,1,1). Greek indices are 
rised and lowere with the Minkowski metric. By convention, Latin indices are raised and lowered with the Kroneker 
symbol that makes no difference between upper and lower Latin indices. 

Bold italic letters (both upper- and lower-case) denote spatial vectors, for instance x = x l = (x l ,x 2 ,x 3 ). Dot and cross 
between two spatial vectors denote the Euclidean scalar and vector products respectively: A B — AjB' = A,B, = A'B 1 , 
and (A x B)' = e ijk AJB k . 

Partial derivative with respect to spatial coordinates x' are denoted as d /dx' or V. Partial derivative with respect to 
four-dimentional coordinates x a are denoted with comma: F a = dF /dx a . Total derivative of function, let say, F(t,x) 
with respect to time is denoted with overdot, that is F = dF/dt, and its partial derivative with respect to time is dF/dt. 
Notice that dF/dt ^ dF/dt. 

G is the universal gravitational constant, and c is the fundamental speed of the Minkowski space-time. 



4. THE IAU 1991 REFERENCE SYSTEMS FRAMEWORK 



Official transition of the astronomical community from the Newtonian to relativistic concepts commenced in 1991 
when a few recommendations (resolution A4) were adopted by the GA of the IAU [20, 97]. 

In the first recommendation, the metric tensor in space-time coordinates (ct,x) centered at the barycenter of an 
ensemble of masses was recommended to be written in the form 
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where c is the fundamental speed of the Minkowski space-time (= speed of light in vacuum), U is the sum of the 
gravitational potentials of the ensemble of masses and of a tidal potential generated by bodies external to the ensemble, 
the latter vanishing at the barycenter. This recommendation recognizes that space-time cannot be described by a single 
coordinate system. The recommended form of the metric tensor can be used not only to describe the barycentric 
celestial reference system (BCRS) of the whole solar system, but also to define the geocentric celestial reference 
system (GCRS) centered on the center of mass of the Earth with a suitable function U, now depending upon geocentric 
coordinates. In analogy to the GCRS a corresponding celestial reference (planetocentric) system can be constructed 
for any other body of the Solar system. 

In the second recommendation, the origin and orientation of the spatial coordinate grids for the BCRS and GCRS 
were defined. The third recommendation defined t = TCB (Barycentric Coordinate Time) and T = TCG (Geocentric 
Coordinate Time) - the time coordinates of the BCRS and GCRS, respectively. The relationship between TCB and 
TCG is given by a full 4-dimensional transformation 
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where x l E = Xe (t ) and v l E are the barycentric coordinate position and velocity of the geocenter, r' E = x' — x' E with x being 
the barycentric position of the observer, and U(t,XE) is the Newtonian potential of all solar system bodies evaluated 
at the geocenter apart from that of the Earth. Fundamental and practical issues associated with the new resolutions on 
the relativistic reference frames and time scales in the solar system were featured in [109-1 12]. 



5. THE IAU 2000 REFERENCE SYSTEMS FRAMEWORK 



In August 2000, new resolutions on reference frames and time scales in the solar system were adopted by the XXIVth 
General Assembly of the IAU [21]. The resolutions are based on the first post-Newtonian approximation of general 
relativity theory and completely abandon the Newtonian paradigm of space and time that changes terminology in 
fundamantal astronomy [1 13]. In this section we reproduce the main results of the BK-DSX theory. For more details, 
the reader is referred to [21] and the original BK-DSX publications. 



5.1. Conventions for the Barycentric Celestial Reference System 

Bary centric Celestial Reference System (BCRS) is denoted x a = (ct,x). BCRS is defined mathematically in terms of 
a metric tensor g a p which components read 
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Here, the post-Newtonian gravitational potential w generalizes the usual Newtonian potential U and w' is the vector 
potential related with gravitomagnetic effects. 

This form of the barycentric metric tensor implies that the barycentric spatial coordinates x' satisfy the harmonic 
gauge condition [92]. The main arguments in favor of the harmonic gauge are: (1) tremendous amount of work on 
isolated astronomical systems in general relativity has been done with the harmonic gauge that was found to be a 
useful and simplifying gauge for various kinds of applications, and (2) in contrast to the gauge adopted in the PPN 
formlism and advocated, for example, by Will in his textbook [43], the harmonic gauge can be defined to higher PN- 
orders by making use of the rather straighforward, post-Newtonian approximations [50], and, in fact, for the exact 
Einstein theory of gravity. 

Assuming space-time to be asymptotically flat (no gravitational fields exist at infinity) in the standard harmonic 
gauge the post-Newtonian field equations of general relativity read 

Dw = -4nGo, (8) 
= -AnG& , (9) 

where □ = -c~ 2 d 2 /dt 2 + V 2 is the (gravity-field) wave operator a = c~ 2 (r 00 + T ss ), a' = c~ l T m , and T^ v are the 
components of the stress-energy tensor of the solar system bodies, T ss = T n + T 22 + r 33 . 

Equations (8), (9) are solved by the post-Newtonian iterations, and in the barycentric coordinate system (t,x l ) one 
has 
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which are to be substituted to the metric tensor (5)-(7). Each of the potentials w and w l can be linearly decomposed in 
two parts 

w = we + w, (12) 
w' = w' E + w' , (13) 

where we and w' E are geopotentials depending on the distribution of Earth's mass and current density expressed in the 
BCRS, and we and w l E are gravitational potentials of all external bodies but the Earth. 



5.2. Conventions for the Geocentric Celestial Reference System 



Geocentric Celestial Reference System (GCRS) is denoted X a = (cT,X). GCRS is defined in terms of the geocentric 
metric tensor G a p with components represented in the following form 
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Here W = W(T,X) is the post-Newtonian gravitational potential in the geocentric system, and W'(T,X) is the 
corresponding vector potential. They satisfy to the same type of the wave equations (8), (9) as the corresponding 
potentials w, and w l in BCRS. 

These geocentric potentials should be split into two parts: potentials We and W E arising from the gravitational action 
of the Earth and external parts w and w l due to dynamic (tidal gravity force) and kinematic (inertial motion) effects. 
The external parts are assumed to vanish at the geocenter and admit an expansion into positive powers of X. Explicitly, 

W(T,X) = W E {T,X)+W^ a (T,X)+W Aya {T,X), (17) 
W\T,X) = W' E {T,X)+WU{T,X)+Wi yn {T,X). (18) 

The Earth's geopotentials W E and W' E are defined in the same way as we and w' E (see equations (10)— (1 1)) but with 
quantities a and a 1 calculated in the GCRS. Wjd n and W^ in are kinematic contributions that are linear in spatial 
coordinates X 

W kia = QiX\ Wi„ = 1 c 2 e ipq (£l p - QP iec ) X" , (19) 

where Qi characterizes the minute deviation of the actual worldline of the origin of the GCRS from geodesic motion 
in the external gravitational field of the solar system bodies 

Q i = d i w(x E )-a E + 0(c- 2 ). (20) 

Here a l E = dv' E /dt is the barycentric acceleration of the origin of the GCRS (geocenter). The function i2p rec describes 
the relativistic precession of dynamically nonrotating spatial axes of GCRS with respect to remote celestial objects 
(quasars) 

^prec^^ 6 ^ [^-\ v id k w{xE)+2d k W } {x E )- l -V j E Q k ^. (21) 

The three terms on the right-hand side of this equation represent the geodetic, Lense-Thirring, and Thomas precessions, 
respectively. One can prove that i2" ner is dominated by geodetic precession amounting to ~ 2" per century plus short- 
periodic terms usually called geodetic nutation. One sees that for Q. 1 = ilp rec the kinematic vector potential W^ in 
vanishes. This implies that dynamical equations of motion of a test body, e.g., a satellite orbiting around the Earth, 
do not contain the Coriolis and centrifugal terms, i.e., the local geocentric spatial coordinates X are dynamically non- 
rotating. For practical reasons, however, the use of kinematically non-rotating geocentric coordinate axes defined by 
condition Q. 1 = is recommended by IAU. Kinematically non-rotating GCRS keeps its spatial axes fixed with respect 
to the axes of the BCRS while dynamically non-rotating GCRS makes spatial axes slowly precessing with respect to 
the BCRS. 

Potentials Wd yn and are generalizations of the Newtonian tidal potential. For example, 

W dyn (T,X) = w(x E +X) -w(x E ) -X i d i w(x E ) + 0(c- 2 ). (22) 

It is easy to check out that a Taylor expansion of w(xe +X) around the point xe gives a polynom starting from the 
quadratic with respect to X terms. We also note that the local gravitational potentials We and W' E of the Earth are 



related to the bary centric gravitational potentials we and w' E by the relativistic transformations [21] 
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where R'j (f ) is a (time-dependent) orthogonal matrix of rotation (with the angular velocity i2p rec ) arising because of the 
choice of the kinematically non-rotating GCRS. In case of the dynamically non-rotating GCRS the matrix R'j(t) = 5j. 



5.3. Transformations between the Reference Systems 

The coordinate transformations between the BCRS and GCRS are found by matching the BCRS and GCRS metric 
tensors in the vicinity of the world line of the Earth by making use of their tensor properties. The transformations are 
written as 
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ci [A{t) + v E -r E ] , c4 



B(t)+B'(t)r' E + B" (t)r E r J E + C(t,x) + O (c' 5 ) , 

0(c- 4 ), 



1 . _ 1 . , 

-v' E v E ■ r E + w{x E )r E + r' E a E ■ r E - -a' E r E 



where r E X x E . and functions A(t),B(t),B'(t),B'j(t),C(t,x) are 
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dt 
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(25) 
(26) 

(27) 
(28) 
(29) 
(30) 
(31) 



Here again x' E , v' E , and a' E are the barycentric position, velocity and acceleration vectors of the Earth, the dot stands 
for the total derivative with respect to t . Let us also remark that the harmonic gauge condition does not fix the function 
C(t,x) uniquely. However, it is reasonable to fix it in the time transformation for practical reasons. 



6. MATCHING THE IAU 2000 FRAMEWORK WITH THE PPN FORMALISM 

General theory of relativity is the most powerful theoretical tool for experimental gravitational physics both in 
the solar system and outside of its boundaries. It passed all tests with unparallel degree of accuracy [114-116]. 
However, alternative theoretical models are required for deeper understanding of the nature of space-time gravitational 
physics and for studying possible violations of general relativistic relationships which may be observed in near-future 
gravitational experiments. 

The goal of this section is to discuss how to incorporate the parametrized post-Newtonian (PPN) formalism [43] 
to the IAU 2000 theory of general relativistic reference frames in the solar system [21]. This will extend domain 
of applicability of the resolutions to more general class of gravity theories. Furthermore, it will make the IAU 2000 
resolutions fully compatible with the JPL equations of motion used for calculation of ephemerides of major planets, 
Sun and Moon. These equations depend on two PPN parameters, j3 and 7 [117] and they are presently compatible 
with the IAU 2000 resolutions only in the case of J3 = 7 = 1 . Rapidly growing precision of optical and radio 
astronomical observations as well as calculation of relativistic equations of motion in gravitational wave astronomy 
urgently demands to work out a PPN theory of relativistic transformations between the local and global coordinate 
systems. 



PPN parameters j3 and y are characteristics of a hypothetical scalar field which perturbs the metric tensor and 
makes it different from general relativity. Scalar fields have not yet been detected but they already play significant role 
in modern physics for many reasons [118, 1 19]. In order to extend the IAU 2000 theory of reference frames to the PPN 
formalism we employed a general class of the scalar-tensor theories of gravity initiated in the pioneering works by 
Jordan [120, 121], Fierz [122] and, especially, Brans and Dicke [123-125]. This class of theories is based on the metric 
tensor g a p representing gravitational field and a scalar field that couples with the metric tensor through the coupling 
function 9(<j>) which we keep arbitrary. Following Will [43], we assume that and 0(0) are analytic functions which 
can be expanded about their cosmological background values and 9. Existence of the scalar field brings about 
dependence of the universal gravitational constant G on the background value of the field which can be considered 
as constant on the time scale much shorter than the Hubble cosmological time. 



6.1. Local Coordinate System in the PPN Formalism 

Our purpose is to develop a theory of relativistic reference frames in an N-body problem (solar system) with two 
parameters j3 and 7 of the PPN formalism. There is a principal difficulty in developing such a theory associated with the 
problem of construction of a local coordinate system in the vicinity of each self-gravitating body (Sun, Earth, planet) 
comprising the N-body system. Standard textbook on the PPN formalism [43] does not contain solution of this problem 
because the original PPN formalism was constructed in a single, asymptotically flat, global PPN coordinates covering 
the entire space-time and having the origin at the barycenter of the solar system. PPN formalism admits existence of 
several fields which are responsible for gravity - scalar, vector, tensor, etc. After imposing boundary conditions on all 
these fields at infinity the standard PPN metric tensor merges their contributions all together in a single expression so 
that they get absorbed to the Newtonian and other general relativistic potentials and their contributions are strongly 
mixed up. Therefore, it is technically impossible in standard PPN formalsim to disentangle the different fields in 
order to find out relativistic space-time transformation between the local frame of a self-gravitating body (Earth) and 
the global PPN coordinates which would be consistent with the law of transformation of the fields imposed by each 
specific theory of gravitation. 

It is quite straightforward to construct the post-Newtonian Fermi coordinates along a world line of a massless 
particle [40]. Such approach can be directly applied in the PPN formalism to construct the Fermi reference frame 
around a world line of, for example, an artificial satellite. However, account for gravitational self-field of the particle 
(extended body) changes physics of the problem and introduces new mathematical aspects to the existing procedure of 
construction of the Fermi frames as well as to the PPN formalism. Only three papers [53-55] have been published 
so far where possible approaches aimed to derive the relativistic transformations between the local (geocentric, 
planetocentric) and the PPN global coordinates were discussed in the framework of the PPN formalism. The approach 
proposed in [53] is based on the formalism that was originally worked out by Ashby and Bertotti [32] and Fukushima 
[42] in order to construct a local inertial frame in the vicinity of a self-gravitating body that is a member of an N- 
body system. In the Ashby-Bertotti formalism the PPN metric tensor is taken in its standard form [43] and it treats all 
massive bodies as point-like monopole massive particles without rotation. Construction of a local inertial frame in the 
vicinity of such massive particle requires to impose some specific restrictions on the world line of the particle. Namely, 
the particle is assumed to be moving along a geodesic defined on the "effective" space-time manifold which is obtained 
by elimination of the body under consideration from the expression of the standard PPN metric tensor. This way of 
introduction of the "effective" manifold is not defined uniquely bringing about an ambiguity in the construction of the 
"effective" manifold that was noticed by Kopeikin [33]. Moreover, the assumption that bodies are point-like and non- 
rotating is insufficient for modern geodesy and relativistic celestial mechanics [62]. For example, planets in the solar 
system and stars in binary systems have appreciable rotational speeds and noticeable higher-order multipole moments. 
Gravitational interaction of the multipole moments of a celestial body with external tidal field does not allow the body 
to move along the geodesic line [33]. Deviation of the body's center-of-mass world line from the geodesic motion 
can be significant and important in numerical calculations of planetary ephemerides [62, 126]) and must be taken into 
account when one constructs a theory of the relativistic reference frames in the N-body system. 

Different approach to the problem of construction of a local (geocentric) reference frame in the PPN formalism was 
proposed in the paper by Klioner and Soffel [54]. These authors have used a phenomenological approach which does 
not assume that the PPN metric tensor in local coordinates must be a solution of the field equations of a specific theory 
of gravity. The intention was to make this kind of formalism as general as possible. To this end, these authors assumed 
that the structure of the metric tensor written down in the local (geocentric) reference frame must have the following 



properties: 

A. gravitational field of external with respect to the Earth bodies (Sun, Moon, planets) is represented in the vicinity 
of the Earth in the form of tidal potentials which should reduce in the Newtonian limit to the Newtonian tidal 
potential, 

B. switching off the tidal potentials must reduce the metric tensor of the local coordinate system to its standard PPN 
form. 

Direct calculations revealed that under assumptions made in [54] the properties (A) and (B) can not be satisfied 
simultaneously. This is a direct consequence of abandoning the field equations for construction of the metric tensor 
in geocentric coordinates and the matching procedure applied in [54] in order to transform the local geocentric 
coordinates to the global barycentric ones. More specifically, at each step of the matching procedure of [54] four 
kinds of different terms in the metric tensors have been singling out and equating separately in the corresponding 
matching equations for the metric tensor: 

• the terms depending on internal potentials of the body under consideration (Earth); 

• the terms which are functions of time only; 

• the terms which are linear functions of the local spatial coordinates; 

• the terms which are quadratic and higher-order polynomials of the local coordinates. 

These matching conditions are implemented in order to solve the matching equations. 

We draw attention of the reader to the problem of choosing the right number of the matching equations. In general 
theory of relativity the only gravitational field variable is the metric tensor. Therefore, it is necessary and sufficient to 
write down the matching equations for the metric tensor only. However, alternative theories of gravity have additional 
fields (scalar, vector, tensor) which contribute to the gravitational field as well. Hence, in these theories one has to use 
matching equations not only for the metric tensor but also for these fields. The paper [54] assumes that it is sufficient 
to solve the matching equations for the metric tensor only in order to obtain complete information about the structure 
of the parametrized post-Newtonian transformation from the local to global frames but this is not correct as shown in 
[55]. In our approach [55] we consistently use the matching equations for both the metric tensor and the scalar field 
which are derived from the field equations. Klioner - Soffel approach to the PPN formalism with local frames taken 
into account admits too many degrees of freedom in matching equations which can not be fixed uniquely. Restriction of 
this freedom can be done ad liberum, thus leading to researcher-dependent ambiguity in the construction and physical 
interpretation of relativistic effects in the local (geocentric) reference frame. 

Our point of view is that in order to eliminate any ambiguities in the construction of the PPN metric tensor in the 
local reference frame of the body under consideration and to apply mathematically rigorous procedure for derivation 
of the relativistic coordinate transformations from the local to global coordinates, a specific theory of gravity must 
be used. This fixes the structure of the field equations and the number of functions entering the PPN metric tensor 
in the local coordinates that becomes exactly equal to the number of matching equations. Hence, all of them can be 
determined uniquely. We propose to build a parametrized theory of relativistic reference frames in the solar system by 
making use of the following procedure: 

1. chose a class of gravitational theories with a well-defined field equations; 

2. impose a specific gauge condition on the metric tensor and other fields to single out a set of global and local 
coordinate systems and to reduce the field equations to a solvable form; 

3. solve the reduced field equations in the global (barycentric) coordinate system x a = (ct,x l ) by imposing fall-off 
boundary conditions at infinity; 

4. solve the reduced field equations in the local (geocentric) coordinate system X a = (cT,X l ) defined in the vicinity 
of the Earth; 

5. make use of the residual gauge freedom to eliminate non-physical degrees of freedom and to find out the most 
general structure of the space-time coordinate transformation between the barycentric and geocentric coordinates; 

6. transform the metric tensor and the other fields from the local coordinates to the global ones by making use of 
the general form of the coordinate transformations found at the previous step; 

7. derive from this transformation a set of matching (the first-order differential and/or algebraic) equations for all 
functions entering the metric tensor and the coordinate transformations; 

8. solve the matching equations and determine all functions entering the matching equations explicitly. 



This procedure works perfectly in the case of general relativity [21] and is valid in the class of the scalar-tensor theories 
of gravity as well [55]. 

The scalar-tensor theory of gravity employed in [55] operates with one tensor, g a p, and one scalar, 0, fields. The 
tensor field g a p is the metric tensor of the Riemannian space-time manifold. The scalar field is generated by matter 
of the bodies comprising an N-body system. We assume that the N-body system consists of extended bodies which 
moves slowly and have weak gravitational field. We use the post-Newtonian approximation (PNA) scheme in order to 
find solutions of the scalar-tensor field equations. It takes into account the post-Newtonian definition of gravitational 
multipole moments introduced by Thorne [127] and further elaborated by Blanchet and Damour [128]. We do not 
specify the internal structure of the bodies so that our consideration is not restricted with the case of a perfect fluid that 
extends the PPN formalism [43]. 



6.2. Transformations between the Reference Systems in the PPN Formalism 

The coordinate transformations between the BCRS and GCRS are found by matching the BCRS and GCRS metric 
tensors and the scalar field in the vicinity of the world line of the Earth. The transformations have the following form 
[55] 
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where = x-X E , and functions A(t),B(t),B'(t),B'j(t),C(t,x) are 
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(36) 
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Here again x' E , v' E , and a' E are the barycentric position, velocity and acceleration vectors of the Earth, the dot stands for 
the total derivative with respect to t . 

As one can see the PPN transformations between the geocentric and barycentric coordinates depend explcitly on 
two PPN parameters j3 and 7. The gauge function C(t,x) remains unaffected by these parameters. In case of general 
relativity, j3 = 7 = 1, and the PPN transformations (32)-(38) are reduced to the form adopted by the International 
Astronomical Union and given in equations (25)— (3 1). 

The PPN geocentric coordinates do not rotate kinematically with respect to the International Celestial Reference 
System (ICRF) defined by a set of primary reference quasars [9, 10]. Dynamically non-rotating geocentric coordinates 
are precessing with respect to ICRF with the angular velocity [55] 



^prec — c 2 



Y+D v e d k w(x E ) + (7+ 1) d k wi{x E ) -\v E Q k 



(39) 



which depends on the PPN parameter 7. Equation (39) is similar to that derived in [92] for a test spinning particle 
(gyroscope) moving in external gravitational field of a rotating mass. We emphasize that equation (39) has been 
derived for a massive self-gravitating body as opposed to the case of the test particle. Equation (39) goes to general 
relativistic expression (21) when 7=1. 



7. MATCHING THE IAU 2000 FRAMEWORK WITH COSMOLOGICAL REFERENCE 

FRAME 



The rapidly growing accuracy of astronomical measurements makes it necessary to take into account some important 
cosmological effects for an adequate interpretation of optical and radio observations of cosmological lenses, anisotropy 
in cosmic microwave background radiation, etc. Some, yet unexplained anomalies in the orbital motion of the solar 
system bodies [23, 24, 129] may be associated with the cosmological expansion and/or other cosmology-related 
effects (multidimensional branes, strings). For this reason, matching of the IAU 2000 framework for reference 
frames in the solar system with the cosmological coordinates of the Friedman-Robertson- Walker (FRW) universe 
[92] becomes vitally important. In this section we outline the main ideas for doing this matching that incorporates 
the Hubble expansion effects to the post-Newtonian approximation of general relativity worked out previously for an 
asymptotically flat space-time. They were worked out in our papers [99, 100]. 

The main observation is that our universe is not asymptotically-flat and, hence, gravitational field does not vanish 
as one approaches infinity. The assumption of vanishing gravitational field and its first derivatives was extensively 
used in the Newtonian and relativistic celestial mechanics of N-body system but, in fact, it is not realistic and 
should be replaced with the cosmological boundary conditions. To this end, we considered a FRW universe driven 
by quintessence (that is, a scalar field imitating the dark matter [119]) 0, and having a spatial curvature equal to zero 
[118, 119]. The background universe is perturbed by a localized, self-gravitating system of N bodies (the solar system) . 
In this model the perturbed metric tensor reads 

gap = a 2 {x){j} a p+h a p) , (40) 

where perturbation h a p of the background FRW metric tensor g a p = a 2 r] a p is caused by the presence of the solar 
system, a{x) is a scale factor of the universe depending on cosmological time T related to coordinate time t by simple 
differential equation dt = a(x)dx. In what follows, a linear combination of the metric perturbations 

yap =h ap_ I 7? «/3/j j (41) 

where h — rj a ^h a p, is more convenient for calculations. 

We discovered a new cosmological gauge [99] which has a number of remarkable properties [130]. In case of the 
background FRW universe with dust equation of state (that is, the background pressure of matter is zero [118, 1 19]) 
this gauge is given by [99, 100] 

f^p=2H(p8 ( f 7 (42) 

where bar denotes a covariant derivative with respect to the background metric g a p, (p = /a 2 , H = d/a is the Hubble 
parameter, and the overdot denotes a time derivative with respect to time x. The gauge (42) generalizes the harmonic 
gauge of asymptotically-flat space-time for the case of the expanding non-fiat background universe. The most attractive 
property of the new gauge is that it allows us to decouple the system of the Einstein equations for the perturbations of 
gravitational field directly in the position space without making their Fourier transform in a series of "plane waves" as 
it is usually done in standard cosmological perturbation theory [119]. 

The gauge (42) drastically simplifies the linearized Einstein equations which are reduced to the following form 

Uy a p-2Hd x y a p+H 2 

where d T = d/dx, □ = -c~ 2 d 2 + V 2 , y= -700 + Ykk, and T a P is the tensor of energy-momentum of matter of the solar 
system. Introducing new notations 700 = Aw/c 2 , yo,- = — Aw'/c 3 , and jij = Aw 1 -' /c 4 , and spiking Einstein's equations in 
time-time, time-space, and space-space components, yield 



U X -2Hd xX + ^H 2 x = -A%Go, (44) 

Dw-2Hd T w = -AnGo-AH 2 x , (45) 

Dw i -2Hd t w i + H 2 w i = —AnGo' , (46) 

Uw ij -2Hd T w ij = -AnGT ij , (47) 



(43) 



where % = w - (p/2, the Hubble parameter H = d/a = 2/t, and densities a = c" 2 (r 00 + T ss ) and & = c^T '. These 
Bessel-wave equations extend the domain of applicability of the post-Newtonian equations (8), (9) to the case of 
expanding universe. 

Equations (44)-(47) consist of eleven equations as besides the metric tensor we have a scalar field that mimics 
a dark matter. First equation (44) describes evolution of the scalar field while the second equation (45) describes 
evolution of the scalar perturbation of the metric tensor. Equation (46) yileds evolution of vector perturbations of the 
metric tensor, and equation (47) describes generation and propagation of gravitational waves by the isolated N-body 
system. Equations (44)-(47) contain all corrections depending on the Hubble parameter and can be solved analytically 
in terms of generalized retarded solution. Exact Green functions for these equations have been found in [99, 100] (see 
also [131]). They revealed that the gravitational perturbations of the isolated system on expanding background depend 
not only on the value of the source taken on the past null cone but also on the value of the gravitational field inside the 
past null cone. Existence of extra terms in the solutions of equations (44)-(47) depending on the Hubble parameter 
brings about cosmological corrections to the Newtonian law of gravity. For example, the post-Newtonian solutions of 
equations (45), (46) with a linear correction due to the Hubble expansion are 

w(r,x) = Gj g( ^ V + §2-^2 J d 3 x'a(r,x')\x-x'\-GH J d\'a{x,x') + 0{ c -") +0(H 2 ) ,(48) 

A*,*) = Gj g '^ V +0(c- 2 )+0(// 2 ). (49) 

Linear with respect to the Hubble parameter, H, term is absent in equation (49) due to the law of conservation of linear 
momentum of an isolated self-gravitating system. Matching these solutions with those defined in the BCRS of the IAU 
2000 framework in equations (10), (1 1) is achieved after expanding all quantities depending on the conformal time x 
in the neighbourhood of the present epoch in powers of the Hubble parameter. 
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